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An  ABCD  beam-propagation  method  was  used  to  build  a  first-order  mathematical  model  of  a  thermal  lens 
effect  from  a  near-infrared  laser  beam  in  water  and  ocular  media.  The  model  was  found  to  fit  experimental 
2-scan  data  best  when  the  thermo-optic  coefficient  Aril  AT  of  liquid  water  at  292  K  was  -4.46  x  10~5  K"1.  The 
physiological  parameters  of  the  human  eye  were  simulated  in  a  simple  eye  model  using  this  fitted  An /AT  value. 
Conservative  model  simulations  for  1150  and  1318  nm  laser  radiation  include  parameter  sets  used  in  experi¬ 
mental  ocular  exposures  performed  by  Zuclich  et  al.  [Health  Phys.  92,  15  (2007)]  to  illustrate  the  transient 
response  of  the  thermal  lens  approaching  the  limits  of  the  retinal  damage  thresholds  for  equivalent  laser  ra¬ 
diation  sources.  ©  2009  Optical  Society  of  America 

OCIS  codes:  330.3350,  190.4710,  120.6810,  170.1020,  170.3660. 


1.  INTRODUCTION 

Laser  safety  for  the  retina  is  established  by  experimental 
studies  at  many  wavelengths  and  exposure  durations  to 
determine  the  smallest  pulse  energy  into  the  eye  that  pro¬ 
duces  a  minimum  visible  lesion  (MVL)  as  seen  with  a  fun¬ 
dus  camera.  Retinal  MVLs  determined  in  the  near- 
infrared  (NIR)  region  are  higher  than  expected  and 
suggest  some  unexplained  laser-tissue  interactions  [1]. 
We  believe  thermal  lensing  plays  a  significant  role  in 
laser-tissue  interactions  influencing  ocular  damage  for 
wavelengths  from  1150-1350  nm,  where  preretinal  me¬ 
dia  strongly  absorb  laser  radiation,  yet  there  is  enough 
energy  remaining  to  pose  a  threat  to  the  sensitive  retinal 
layers. 

Thermal  retinal  damage  models  use  laser-induced  heat 
generation,  heat  conduction,  and  the  Arrhenius  damage 
integral  to  predict  the  formation  of  a  laser-induced  lesion 
in  the  retina  [2] .  These  models  require  an  understanding 
of  laser  beam  spot  size  and  energy  at  the  retina  that  de¬ 
pend  on  optical  properties  of  the  preretinal  media  includ¬ 
ing  dispersion,  absorption,  and  physical  dimensions  of  the 
eye.  Similar  fundamental  understandings  are  necessary 
to  estimate  an  action  spectra  used  to  make  predictions  on 
general  trends  by  fitting  available  MVL  data  over  the  en¬ 
tire  UV-visible-IR  spectrum  [1,3—5].  To  date,  thermal 
lensing  has  been  a  relatively  unexplored  laser-tissue  in¬ 
teraction  in  the  eye  and  has  not  been  incorporated  in  any 
eye  model  used  to  predict  retinal  damage  thresholds. 

The  discovery  of  thermal  lensing  by  Gordon  et  al.  in 
1965  led  to  the  development  of  numerous  mathematical 
models  [6].  The  phenomenon  is  the  result  of  a  tempera¬ 
ture  gradient,  typically  assumed  to  be  radially  symmetric, 
formed  by  the  absorption  of  laser  light  in  the  medium.  As 
the  temperature  T  of  the  medium  increases,  the  local  den- 
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sity  p  decreases.  This  leads  to  a  decrease  in  the  index  of 
refraction  n  resulting  in  the  formation  of  a  negative  lens. 
The  temperature  gradient  is  shaped  by  the  beam  profile 
and  the  thermal  diffusivity  of  the  medium  [7-10].  The 
creation  of  a  thermal  lens  in  ocular  media  means  the  spot 
size  formed  at  the  retina  changes  dynamically  as  a  func¬ 
tion  of  the  coupled  transient  response  of  heat  generated 
by  the  beam  profile  and  thermal  diffusion. 

Many  thermal  lensing  models  begin  with  a  solution  of 
the  heat  conduction  equation  in  combination  with  an 
ABCD  ray-tracing  method  [9—11] .  The  model  assumes  a 
radially  symmetric  temperature  gradient  in  an  absorbing 
medium  where  no  heat  is  exchanged  along  the  optical  axis 
z. 

A  closed-aperture  z-scan  is  one  technique  to  correlate 
transmission  data  through  an  absorbing  medium  to  a 
thermal  lens  gradient  [12].  A  single-beam,  closed- 
aperture  z-scan  using  1313  nm  light  provides  an  estimate 
of  thermal  lensing  in  water.  The  thermal  characteristics 
of  water  are  well  documented  and  show  similar  linear  ab¬ 
sorption  coefficients  to  those  of  ocular  media  at  NIR  wave¬ 
lengths,  making  the  z-scan  a  basic  method  for  initial 
evaluation  of  the  thermal  lensing  model  [13-17], 

In  this  paper,  data  from  a  single-beam,  closed-aperture 
z-scan  using  a  continuous  wave  (CW)  1313  nm  laser  ra¬ 
diation  source  are  presented  and  compared  with  a  first- 
order  model  for  thermal  lensing.  Data  from  the  z-scan  are 
used  to  determine  an  estimate  of  the  change  in  index  of 
refraction  with  the  change  in  temperature  dn/dT  for  this 
region  of  laser  radiation.  Having  demonstrated  a  reason¬ 
able  fit  to  z-scan  data,  the  first-order  model  is  then  used 
as  a  tool  for  predicting  trends  of  the  thermal  lensing  effect 
in  a  virtual  eye.  Preliminary  data  relating  to  this  work 
have  been  reported  [18-21],  The  first-order  model  pre- 
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sented  here  does  not  include  heat  conduction  along  the 
z-axis  nor  the  effects  of  convection  in  the  absorbing  me¬ 
dia.  The  more  viscous  properties  of  ocular  media  com¬ 
pared  with  water  preclude  significant  convectional  flow. 


2.  BACKGROUND 

A.  Z-Scan  Spectroscopy 

The  z-scan  is  a  common  analytical  technique  used  to  mea¬ 
sure  the  thermal  lens  effect.  The  technique  is  a  well  es¬ 
tablished  method  in  chemistry  for  chromatography  and 
electrophoresis  analyses  and  in  electro-optics  for  testing 
the  characteristics  of  optical  limiters  [7—12]. 

A  single-beam  z-scan  is  the  simplest  of  z-scan  tech¬ 
niques,  where  a  single  laser  is  used  as  the  pump  and 
probe.  The  laser  beam  is  collimated  and  aligned  to  be  sent 
through  a  focusing  lens.  A  chopper  or  shutter  placed  after 
the  focusing  lens  is  traditionally  used  as  the  control  for 
the  exposure  duration  [12].  The  sample  is  placed  at  posi¬ 
tion  z  in  the  beam’s  path  after  the  shutter  and  at  an  angle 
just  a  few  degrees  from  the  z  axis  to  prevent  reflective  in¬ 
terference  from  the  cuvette  windows.  A  detector  is  placed 
in  the  far  field  from  the  focusing  lens.  An  aperture  is 
placed  just  in  front  of  the  detector,  giving  the  nomencla¬ 
ture  of  a  “closed-aperture  z-scan.”  Irradiance  data  are  col¬ 
lected  as  a  function  of  time  as  the  sample  is  held  in  a  sta¬ 
tionary  position.  The  sample  is  allowed  time  to  relax  back 
to  its  initial  thermal  state,  then  moved  to  a  new  z  position 
before  beginning  the  next  exposure.  Conventional  z-scan 
data  report  the  z-position  of  the  sample  in  relation  to  the 
distance  from  the  focal  plane  of  the  focusing  lens.  The  ab¬ 
sorption  of  laser  radiation  in  liquid  media  causes  the  for¬ 
mation  of  a  negative  thermal  lens.  The  irradiance  mea¬ 
surements  are  analyzed  at  a  selected  time  and 
normalized  to  a  reference  irradiance  signal  taken  for  the 
sample  placed  at  a  far  distance — approximately  one  focal 
length — from  the  focus  of  the  lens  (Fig.  1).  The  z-scan 
data  [Fig.  1(d)]  depict  an  increase  [Fig.  1(b)]  and  decrease 
[Fig.  1(c)]  in  relative  irradiance  for  sample  positions  be¬ 
fore  and  after  the  focus,  respectively,  with  an  inflection 
point  occurring  at  the  focal  plane  of  the  lens  [Fig.  1(a)]. 

Selections  of  the  lens  and  the  sample  thickness  are  im¬ 
portant  for  the  z-scan.  The  irradiance  of  the  thermal  lens 
is  largely  related  to  the  generation  of  heat  [W/m3]  in  the 
sample,  making  the  choice  of  a  short  focal  length  lens 
seem  desirable.  However,  since  the  thermal  lens  measure¬ 
ment  is  very  sensitive  to  confocal  distance,  Franko  and 
Tran  describe  using  a  sample  thickness  Ts  that  is  at  least 
one-third  or  less  than  the  confocal  distance: 


127 T(0„ 


To  — - 


(1) 


where  u>min  is  the  minimum  1/e2  beam  waist  radius  of  the 
probe  beam  with  wavelength  X  [12].  The  minimum  waist 
can  be  approximated  as  fk/2d  where  d  is  the  beam’s  1/e2 
diameter  at  the  focusing  lens  with  focal  length  f.  As  a  re¬ 
sult,  the  selected  sample  path  length  is  proportional  to 
the  confocal  distance,  or  twice  the  Rayleigh  range,  of  the 
system  design. 

A  single-beam,  closed-aperture  z-scan  allows  for  the 
transient  information  of  only  the  excitation  of  the  sample 


(a) 


Focusing  laser  beam 
z  =■  0  at  focus 
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Detector 

/ 
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- 1 - ► 

0 
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Fig.  1.  Example  of  a  negative  thermal  lens  z-scan.  The  cases  for 
the  samples  exposed  at  the  z  positions:  (a)  very  far  from  the  focus 
of  the  lens  and  very  near  or  at  the  focus  resulting  in  unity  nor¬ 
malized  irradiance,  (b)  before  the  focus  resulting  in  an  increase 
in  normalized  signal  irradiance,  and  (c)  after  the  focus  resulting 
in  a  decrease  in  signal  irradiance.  (d)  Resulting  z-scan  data  set. 


[12],  When  both  the  excitation  and  relaxation  of  the 
sample  are  desired,  a  dual-beam  z-scan  should  be  used. 


B.  Mathematical  Model 

Using  Yariv’s  approach  for  deriving  the  relationship  of  a 
medium  with  a  quadratic  index  of  refraction  to  a  ray,  the 
index  of  refraction  n  can  be  described  as  a  function  of  ra¬ 
dial  distance  r  from  the  optical  axis  z  and  temperature  T 
[11],  Using  a  Maclaurin  series  to  expand  an  expression  for 
n{r,T)  gives 


j  cfn  \  r2  (  cPn 

♦iiU* 


(2) 


r=0 


Truncating  the  series  and  applying  the  chain  rule  gives 
(  dn  cfT\ 


n(r,T)  =  n0  +  r  - 


-=o 


r2  d  /  dn  dT 
2  Or  \rfT  <tr 


(3) 


r=0 


The  change  in  temperature  with  respect  to  radial  distance 
at  radial  position  r=0  can  be  set  to  zero,  simplifying  the 
expression  to  just  the  background  index  of  refraction  n0 
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and  the  last  term  containing  the  square  of  the  radius  r. 
From  this  point,  applying  the  product  rule  gives  an  ex¬ 
pression  that  is  then  evaluated  for  r/T/rfr-*0  at  r=0, 
where  the  maximum  irradiance  of  the  beam  is  found,  re¬ 
vealing  the  expression 


n{r  ,T)  =n0  +  - 


d  dn 


ar\ar 


ar 


drr~°  2 


a  /  aT 


t/r\  dr 


dn 


ar 


(4) 


which  further  reduces  to 


r2  cPT  dn 


(5) 


The  background  index  of  refraction  n0  must  be  selected 
for  the  corresponding  medium  and  wavelength  of  interest. 
The  term  iPTh'tr 2  is  solved  using  a  Green’s  function  solu¬ 
tion  to  the  heat  conduction  equation,  explained  below 
[8-10],  For  cases  where  the  laser  pulse  is  longer  than  the 
acoustic  transit  time,  as  is  the  case  with  CW  laser  expo¬ 
sures,  the  change  in  n,  An,  can  be  approximated  by 
(dnlAT)(AT)  [7],  The  term  dn/d7’  is  referred  to  as  the 
thermo-optic  coefficient  and  is  wavelength  dependent 
[8,22]. 

An  equivalent  expression  for  n(r,D  for  a  medium  with 
a  quadratic  index  of  refraction  is  described  by  Yariv  as 
[11] 

n  =  n0(  1  -  —  I  •  (6) 


A  solution  for  X  is  determined  by  considering  only  CW 
cases  and  setting  Yariv’s  equation  for  n,  Eq.  (6),  equal  to 
the  solution  for  n(r,T),  Eq.  (5),  giving 


1  dn  rfT 
n0  AT  dr 2 


(7) 


The  term  X  is  needed  to  describe  how  a  ray  would  move 
through  a  quadratic-index  medium  using  an  ABCD 
propagation  matrix.  The  ABCD  matrix  for  propagating  a 
beam  through  a  negative-lenslike  medium  is  given  as  [11] 


A 

B~ 

cosh(Y1/2Az) 

XV2  sinh(ZK2A z) 

C 

D 

Xm  sinh(X1/2Az) 

cosh(Y1/2Az) 

The  term  A z  is  the  axial  step  or  length  of  the  step  through 
which  the  beam  is  propagating.  A  complex  beam  param¬ 
eter  q(z )  is  used  to  describe  the  ray  at  position  z  anywhere 
along  the  z  axis  in  a  medium  as 


Aq(z)0id  +  B 
Cq(z)old  +  D' 


(9) 


The  complex  beam  parameter  is  described  by  the  real  and 
imaginary  portions  of  a  wavefront  by 


versing  the  medium  and  the  1/e2  beam  radius  co(z).  The 
value  of  m2,  or  m -squared,  describes  the  quality  of  colli- 
mation  of  the  beam  [23,24],  Ideally,  a  Gaussian  beam 
would  have  an  m2  value  of  1.0,  but  realistically,  the  value 
is  greater  than  1.0. 

To  solve  for  the  cPT/ dr 2,  the  eye  can  be  treated  as  a  cy- 
lindrically  symmetric  medium  with  thickness  z.  For  our 
first-order  model,  thermal  diffusion  along  the  optical  axis 
is  assumed  to  be  zero  (<?T/ofe= 0).  Thus,  we  assume  diffu¬ 
sion  occurs  along  the  radial  direction.  However,  the  heat- 
source  term  S(z)  is  attenuated  in  the  z  direction  according 
to  Beer’s  Law  [see  Eq.  (15)  below].  The  thermal  diffusion 
equation, 


l  a 

r  ifr 


aT 

+  S(z)  =  pcp-—. 


(ID 


must  be  solved  first  by  differentiating  the  heat  equation 
with  respect  to  T(r,t),  then  integrating  to  find  < PT/dr2. 
Here  t  is  time,  r  is  the  magnitude  of  the  radial  vector  de¬ 
scribed  by  position  (x,y),  t)  is  the  thermal  diffusivity  in 
cm2/s,  and  S(z)  is  the  source  term.  Note:  y=K/pcp  where  k 
is  the  thermal  conductivity  in  W/cm  K,  p  is  the  density  in 
g/cm3,  and  cp  is  the  specific  heat  at  constant  pressure  in 
J/gK. 

We  modeled  our  eye  as  an  infinite  cylinder  with  heat 
generation  originating  from  the  Gaussian  beam  launched 
into  the  eye.  Since  the  conduction  of  heat  is  assumed  to  be 
purely  a  radial  process,  the  diffusion  equation  is  indepen¬ 
dent  of  angle  and  axial  position,  allowing  the  radial  vector 
r  to  reduce  to  a  scalar  radial  magnitude  r  [25,26].  Apply¬ 
ing  the  product  rule  and  rearranging  variables  in  Eq.  ( 11) 
gives  the  heat  equation 


ipT  1  aT  1 


1  ar 


__  +  —  +  -S(z)  -  — . 

■q  M 


dr i 


r  dr  k 


(12) 


The  Green’s  function  for  a  virtual  cylindrical  eye  with 
no  diffusion  along  z  is  taken  as 


G(r,t;r',t')  =  - 


4Trr/(t  - 1') 

(r2  +  r'2) 


Xexp 


4  y(t-t') 


2v(t-t') 


,  (13) 


where  r'  describes  the  position  of  the  cylindrical  surface 
heat  source,  and  70  is  the  modified  Bessel  function  of  the 
first  kind  [10,26].  The  units  of  this  Green’s  function  are  in 
cm-2.  G(r,t;r',t')  = 0  for  all  positions  where  r  is  on  the 
surface  of  the  cylindrical  medium.  Since  only  CW  cases 
are  being  considered  here,  f'  =  0.  Applying  the  Green’s 
function  Eq.  (13)  to  the  heat  diffusion  Eq.  (12),  for  a  sym¬ 
metrical  Gaussian  beam  results  in 


1  1  m2\ 

q{z)  R(z)  '  mo(z)2’  ^  ^ 

where  the  real  portion  of  the  complex  beam  parameter 
contains  information  about  the  radius  of  curvature  R(z), 
while  the  imaginary  portion  contains  the  wavelength  tra- 


T(z,r,t) 


1 

pcP 


r'  dr'  S(z,r')G(r,t;r') . 


(14) 


The  expression 
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S(z,r') 


2  HaPi 

TTCO2 


exp 


-  2(r')2 


(15) 


describes  the  source  term  S(z,r')  for  the  symmetrical 
Gaussian  beam;  for  the  eye  in  W/cm3  it  can  be  found  by 
applying  Beer’s  law  assuming  a  Gaussian  beam  [9,10],  In 
Eq.  (15),  fia  is  the  linear  absorption  coefficient  in  cm-1,  u> 
is  the  1/e2  beam  radius,  and  r'  denotes  the  location  of  the 
heat  source.  Pz  is  the  power  in  watts  delivered  at  position 
z  in  the  sample,  given  by  Beer’s  law  as  Pz=Pp  exp (-/xaz), 
with  P0  as  the  initial  power  delivered  at  the  front 
surface  of  the  sample.  The  power  is  expressed  in  an 
array  of  values  by  Beer’s  law  through  P2+ii2 
=P0  exp(-naz)exp(-/j.aAz),  where  z  is  the  position  along 
the  z  axis  within  the  sample  and  A z  is  the  same  spatial 
step  described  earlier  in  Eq.  (8). 

Using  Mathematics  5.2,  a  solution  for  iPTI dr2  was  veri¬ 
fied  to  be 


<PT 

fir2 


-  v8/J.lrP2 
2 

7TKC0 


t 

8r/t  +  (x>2 


(16) 


The  beam  and  subsequent  radial  temperature  profiles  are 
assumed  to  be  Gaussian  over  all  time  t.  The  critical  time 
is  defined  by  tc=u>2/8r/,  allowing  the  solution  to  reduce  to 

-n*Pz 

fir2  7 tio2k 
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(17) 


Analysis  of  the  problem  can  include  observing  the  tem¬ 
perature  rise  along  the  optical  axis.  Given  our  assumption 
of  fiT/fiz=0  and  initial  condition  fiT/fit= 0,  the  tempera¬ 
ture  rise  along  the  z  axis  is  representative  of  the  central 
temperature  resulting  from  diffusion  along  the  radius  r. 
The  temperature  rise  along  the  z  axis  at  r  =  0  is  then  de¬ 
termined  by  evaluating  the  heat  conduction  equation  at 
r=0,  giving  [this  requires  applying  L’Hopital’s  Rule  to  the 
term  (l/r)d77dr] 


=  V 


r= 0 
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‘  fir2 


- S(z ) 


(18) 


The  source  term  at  r=0  simplifies  to  S(z)  =  2/j.aPz/wu)2, 
where  Pz  was  previously  described.  Numerical  integration 
of  Eq.  ( 18)  with  respect  to  time  can  then  be  used  to  deter¬ 
mine  an  approximate  solution  to  the  temperature  T.  We 
assume  a  uniform  temperature  T0  when  /0=0  with  an  ini¬ 
tial  beam  radius  io0  created  from  the  optical  geometry  of 
the  system  with  no  thermal  lensing.  Methods  for  solving 
this  integral  have  been  reported  in  Swofford  and  Morrell 
[27]. 


3.  MATERIALS  AND  METHODS 

A.  Single-Beam,  Closed-Aperture  Z-Scan 
A  1313  nm  IRCL-150-1313-P-L  CrystaLaser  was  colli¬ 
mated  and  aligned  (Fig.  2)  to  create  a  1/e2  beam  diameter 
at  the  focusing  lens  of  3.5  mm  and  an  m -squared  of  1.16. 
A  400  mm  BK7  planoconvex  lens  was  used  as  the  focusing 
lens.  A  knife-edge  test  found  the  1313  nm  laser  focus  to  be 
42  cm  from  the  center  of  the  lens. 

The  power  of  the  1313  nm  laser  was  attenuated  by  us¬ 
ing  a  Thorlab  NEK01  neutral  density  (ND)  filter  set.  A 


i2  li 
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1  L3  s 
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Fig.  2.  Closed-aperture,  single-beam  z-scan  set  up.  The  z  axis  is 
referenced  to  the  focal  position  of  the  lens  found  using  a  knife- 
edge  technique  in  air,  with  no  cuvette  in  the  beam  path.  The  la¬ 
ser  was  a  1313  nm  CrystaLaser  IRCL-150-1313-P-L.  M1-M4,  IR- 
coated  mirrors;  LI  and  L2,  IR-coated  lenses  created  a  beam 
expander  of  =1.5X;  L3,  400  mm  ThorLabs  planoconvex  1  in. 
lens;  SW,  1  cm  diameter  Uniblitz  shutter;  Quartz  cuvette  from 
Stama  Cells;  A,  Aperture  set  to  2.5  mm  diameter;  D,  818-IR 
Newport  detector. 


Uniblitz  10  mm  shutter  was  used  to  control  the  delivery 
of  the  laser  beam  to  the  sample.  The  focusing  lens  and 
any  ND  filters  were  placed  before  the  Uniblitz  shutter  to 
allow  for  any  thermal  lens  in  the  optical  components  to 
remain  in  steady  state  throughout  the  duration  of  the  ex¬ 
periment. 

All  experiments  were  run  at  room  temperatures  be¬ 
tween  17-20  °C.  A  sample  of  ultrapure  deionized  water 
was  placed  in  a  10  mm  quartz  cuvette  purchased  from 
Starna  Cells  Inc.  Samples  were  run  in  sets  at  100  and 
1000  ms  exposures  for  powers  of  2.60  and  48  mW.  The 
linear  absorption  coefficient  of  water  for  1313  nm  was  es¬ 
timated  to  be  1.27  cm"1  [13], 

A  Newport  818-IR  detector  was  placed  117  cm  down  the 
optical  path  from  the  focusing  lens.  A  2.5  mm  aperture 
placed  in  front  of  the  detector  was  found  to  yield  the  larg¬ 
est  dynamic  range  in  normalized  irradiance  data  for  our 
z-scan  system.  A  detailed  analysis  of  the  selection  of  an 
aperture  has  been  reported  by  Sheik-Bahae  et  al.  [28]. 
The  irradiance  signal  from  the  detector  was  sent  to  a 
Melles  Griot  large  gain  transimpedance  amplifier  before 
being  recorded  in  Labview,  A  Stanford  Research  Systems 
signal  delay  generator  was  used  to  synchronize  the  data 
acquisition  with  the  shutter.  The  sample  was  then  placed 
on  a  stage  along  the  optical  path  of  the  focusing  lens  rang¬ 
ing  from  -18  to  22  cm  in  1  cm  increments  relative  to  the 
focal  position  of  the  lens.  The  irradiance  values  were  nor¬ 
malized  by  values  measured  for  a  sample  position  at  least 
22  cm  from  the  focus. 

In  earlier  reports  of  our  first-order  thermal  lensing 
model,  we  used  constant  thermal  coefficients  and  discov¬ 
ered  a  strong  sensitivity  in  our  early  model  to  small 
changes  in  these  coefficients  [18],  To  model  the  z-scan  ex¬ 
periments  presented  here,  values  for  thermal  conductiv¬ 
ity  k,  specific  heat  at  constant  pressure  cp,  and  density  p 
as  a  function  of  temperature  T  were  obtained  from  Wag¬ 
ner  and  Kruse  at  0.1  MPa  [15].  Data  between  10-90  ”C 
were  fitted  to  a  polynomial  in  excel  to  obtain  continuous 
temperature-dependent  thermal  properties.  These  poly¬ 
nomials  are  given  in  Table  1.  Thermal  diffusivity  rj  was 
calculated  from  r)=xl pcp,  then  fitted  to  a  polynomial.  The 
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Table  1.  Thermal  Conductivity  k,  Density  p,  Spe¬ 
cific  Heat  at  Constant  Pressure  cp,  and  Thermal 
Diffusivity  r/,  As  a  Function  of  Temperature  T 
Found  to  Fit  Data  for  Liquid  Waterafor  0.1  MPa 
between  283.15  and  363.15  K 


Coefficient  As  a  Function  of 

Symbol  [Units] 

Temperature,  T  [K] 

k  [W/(cm  K)]  = 

-9.99  x  10“3T2  +  7.649  x  l(r6T. . . 

-7.851  xio-3 

p  [g/cm3]  = 

-3.461  x  10'6T2+ 1.798  x  10-3T. . . 
+7.684  xlO-1 

1.743  xlO-9T4- 2.322  Xl0-6103... 

Cp  [J/(g  K)J  = 

+  1.169  X10-3T2-2.632X10-1T... 

+2.649  X101 

t]  [cm2/s]  = 

-2.207  X  10-3T2  + 1.778  x  10-5T. . . 

-1.882  xlO-3 

"From  Wagner  and  Kruse  [15]. 


linear  regression  R2  value  was  greater  than  0.995  for  all 
polynomial  fits  to  the  Wagner/Kruse  data. 

B.  Model  of  the  Virtual  Eye 

The  eye  is  a  complex  optical  system,  but  for  our  purposes 
in  assessing  trends  of  thermal  leasing  in  an  eye,  we  have 
simplified  our  model  eye  down  to  a  basic  compartmental 
system  consisting  of  a  single  curved  surface.  Limitations 
regarding  the  simple  eye  model  are  reserved  for  the  dis¬ 
cussion.  Thermal  lensing  in  a  simple  eye  model  was  simu¬ 
lated  for  two  wavelengths:  1318  nm,  as  this  was  the 
wavelength  used  in  retinal  ED50  studies  by  Zuclich  et  al. 
[3],  and  1150  nm,  to  contrast  variations  in  thermal  lens¬ 
ing  from  a  change  in  absorption  coefficients.  The  eye  was 
modeled  as  a  four-chamber  compartment  partitioned  into 
cornea,  aqueous,  lens,  and  vitreous  chambers  with  a  ra¬ 
dius  of  curvature  existing  only  on  the  front  surface  of  the 
air-cornea  boundary.  This  radius  of  curvature  is  essential 
in  the  ABCD  beam  propagation  for  our  thermal  lensing 
model.  Linear  absorption  coefficients  for  each  of  these 
ocular  media  for  1150  and  1318  nm  are  in  Table  2.  A  re¬ 
view  of  the  literature  found  the  average  thickness  of  these 
chambers  for  a  human  to  be  0.05,  0.31,  0.36,  and  1.72  cm 
for  cornea,  aqueous,  lens  and  vitreous,  respectively  (Fig. 
3)  [29-31].  With  a  known  refraction  and  physiological  dis¬ 
tance  to  a  virtual  retinal  plane  of  2.44  cm,  a  radius  of  cur¬ 
vature  rc  of  0.61  cm  was  found  to  set  the  589  nm  wave¬ 
length  to  a  minimum  focus  at  z  =  2.44  cm. 

Refractive  error  measurements  for  wavelengths  beyond 
900  nm  in  the  eye  are  limited,  thus  chromatic  dispersion 
for  wavelengths  greater  than  900  nm  were  based  on  data 
trends  in  Fernandez  et  al.  [32].  Using  these  trends  in  re¬ 
fractive  error,  a  bulk  refractive  index  value  as  a  function 


Virtual 

Retina 


Fig.  3.  Simple  eye  model  with  a  radius  of  curvature  of  0.61  cm 
on  the  front  surface  of  the  cornea  brings  a  collimated  589  nm 
beam  to  a  minimum  at  the  virtual  retinal  plane  2.44  cm  after  the 
front  surface  of  the  eye.  The  vitreous  layer  is  infinitely  long,  with 
beam  trend  observations  made  at  the  position  of  the  virtual 
retina. 

of  wavelength  was  determined  by  a  Sellmeier  fit,  as  de¬ 
scribed  in  an  earlier  publication  by  Vincelette  et  al.,  to  ac¬ 
count  for  chromatic  dispersion  in  our  eye  model  [1],  For  a 
given  wavelength,  the  initial  values  of  u>(z)  in  the  virtual 
eye  are  calculated  from  the  ABCD  propagation  method 
described  by  Yariv  from  Eqs.  (6M10)  using  the  corre¬ 
sponding  bulk  refractive  index  from  the  Sellmeier  fit  re¬ 
ported  in  Vincelette  et  al.  [1].  The  thermal  conductivity  k 
specific  heat  at  constant  pressure  cp,  and  density  p  of  ocu¬ 
lar  media  were  held  as  constants  reported  by  Okuno  et  al. 
[33], 

C.  First-Order  Thermal  Lensing  Mathematical  Model 
The  mathematical  model  presented  earlier  was  built  in 
Matlab  to  predict  the  thermal  lensing  effect  in  a  single¬ 
beam,  closed-aperture  z-scan  and  in  the  reduced  eye 
model.  The  thermal  lensing  model  requires  an  initial  seed 
of  1/e2  beam-waist  values  based  on  no  thermal  lensing; 
these  become  the  beam-waist  values  at  t =0  s.  Subsequent 
1/e2  waist  values  for  f>0  are  then  found  through  an  it¬ 
erative  process  using  small  time  steps  initially,  gradually 
increasing  the  time  step  size  as  the  thermal  lens  reaches 
steady-state. 

A  value  for  dn/dT  was  selected  for  the  model  based  on 
a  best  fit  to  the  z-scan  data.  This  value  of  dn/dT  was  then 
used  in  modeling  the  thermal  lensing  effect  in  the  simple 
eye  model. 

To  determine  the  thermo-optic  coefficient  from  litera¬ 
ture,  the  index  of  refraction  data  at  0.1  MPa  (~1  atm)  for 
laser  radiation  at  404,  589,  632, 1013  and  2325  nm  in  wa¬ 
ter  in  the  liquid  phase  between  10  and  90  °  C  as  reported 
by  Schiebener  et  al.  were  fitted  to  third-order  polynomials 
in  EXCEL  to  determine  tt(T)x  [16].  The  first  derivative  was 


Table  2.  Linear  Absorption  Coefficients  for  Ocular  Media  and  Water  at  1150  and  1318  nm“ 


Absorption  Coefficient,  pa  [cm-1] 

[nm] 

Cornea 

Aqueous 

Lens 

Vitreous 

Water 

1150 

1.51 

1.22 

0.61 

1.13 

1.78 

1318 

1.98 

2.01 

0.97 

1.90 

1.34 

'From  [13,17], 
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Fig.  4.  Values  of  dn/dT  (in  K-1)  as  a  function  of  temperature  for 
water  in  the  liquid  phase  for  five  different  wavelengths.  Values 
here  were  based  on  taking  the  derivative  of  n{T)P  at  each  of  the 
five  wavelengths  as  reported  by  Schiebener  et  al.  at  0.1  MPa  [16], 

then  taken  from  the  fit  to  determine  dn/dT(T)  for  five  dif¬ 
ferent  wavelengths  as  shown  in  Fig.  4. 

4.  RESULTS 

A.  Model  Compared  with  Experimental  Z-Scan  Data 
By  substituting  Eq.  (17)  into  Eq.  (7)  and  using  the  ABCD 
beam-propagation  method  described  in  Eqs.  (8>— (10)  the 
model  was  arranged  to  simulate  the  2-scan  set  up  de¬ 
picted  in  Fig.  2.  A  comparison  of  the  results  from  this  por¬ 
tion  of  the  model  to  experimental  2 -scan  data  for  varying 
input  power  and  a  10  mm  thick  water-filled  cuvette  at 
room  temperature  can  be  found  in  Figs.  5  and  6.  The 
2-scan  data  presented  here  reflect  the  average  of  three  ex¬ 
perimental  runs.  The  standard  deviation  of  the  normal¬ 
ized  intensity  (for  experimental  data  results)  was  ap¬ 
proximately  ±0.02  arbitrary  units. 

A  bulk  value  of  -4.46  X  10~5  K-1  for  dn/dT  was  found  to 
give  the  best  fit  to  the  experimental  1313  ran  radiation 
2-scan  data  in  terms  of  line  width  before  and  after  the  fo¬ 
cal  position  at  2=0  and  yielded  the  best  approximation  to 


the  slope  between  the  maxima  and  minima.  A  detailed 
discussion  of  the  impact  of  the  value  of  dn/dT  is  reserved 
for  later. 

B.  Thermal  Lensing  Model  Trends  in  Virtual  Eye 
The  dn/dT  value  found  to  fit  the  2-scan  data  was  used  in 
the  simple  eye  model  presented  in  Fig.  3  to  make  predic¬ 
tions  on  the  trends  of  thermal  lensing  in  a  human  eye. 
Analysis  of  the  Schiebener  et  al.  data  [16]  (Fig.  4)  sug¬ 
gests  that  this  value  of  dn/dT  will  approximately  double 
when  the  ambient  temperature  is  increased  from  room  to 
body  temperature,  lb  test  for  the  significance  in  this 
variation  due  to  an  increasing  thermo-optic  coefficient  in 
our  model,  we  used  a  constant  dn/dT  value  of  -4.46 
X  10'5  K_1  for  the  results  in  Figs.  7-10  and  then  doubled 
dn/dT  to  -8.92  x  10~5  K~x  for  the  results  in  Fig.  11.  Again, 
Eq.  (17)  was  substituted  into  Eq.  (7)  and  the  ABCD  beam- 
propagation  method  described  in  Eqs.  (8)-(10)  was  used, 
except  now  the  model  was  arranged  to  simulate  the  vir¬ 
tual  eye  (Fig.  3).  All  other  values  for  the  simple  eye  model 
have  been  previously  described. 

The  influence  of  power  from  a  Gaussian  1318  nm  laser 
radiation  source  input  into  the  eye  on  the  beam’s  radius 
at  the  retina  as  a  function  of  time  is  shown  in  Fig.  7.  The 
results  in  Figs.  7  and  9-11  were  normalized  to  the  initial 
(  =  0  beam  radius  at  the  retina  in  order  to  examine  the 
magnitude  of  the  increase  in  the  beam’s  radius  (i.e.,  a  fac¬ 
tor  of  two,  ten,  etc.)  due  to  thermal  lensing.  Where  appro¬ 
priate,  the  value  for  normalizing  the  results  is  given  in 
the  caption.  The  20  W  case,  shown  in  Fig.  7,  represents 
the  upper  limit  of  the  1318  nm,  CW  laser  radiation  source 
that  was  used  in  CW  nonhuman  primate  retinal  damage 
studies  conducted  by  Zuclich  et  al.  [34] . 

Keeping  the  power  delivered  to  the  cornea  at  3  W,  the 
influence  of  the  beam  diameter  from  a  Gaussian  1318  nm 
laser  radiation  source  input  at  the  cornea  on  the  beam  ra¬ 
dius  at  the  retina  as  a  function  of  time  is  shown  in  Fig.  8. 
Parameters  were  for  a  3  W,  CW,  1318  nm  laser  radiation 
source  delivering  collimated  1/e2  Gaussian  beam  radii  of 
0.5,  0.75,  1.25,  1.75,  and  2.25  mm  at  the  virtual  cornea. 
Note  that  Zuclich  et  al.  reported  the  corneal  ED50  for  a 
CW,  1318  nm  laser  radiation  delivered  to  the  cornea  for 
0.28  s  exposures  is  0.56  J  [34],  Zuclich  et  al.  used  a  beam 


(a)  (b) 


Fig.  5.  2-scan  data  for  a  10  mm  water-filled  cuvette  exposed  to  1313  nm  CW  laser  source  for  1  s  at  (a)  2.6  mW  and  (b)  48  mW.  The  time 
response  for  the  maxima  and  minima  from  each  peak  can  be  found  in  Fig.  6. 
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Time  (s)  Time  (s) 


Fig.  6.  Transient  response  for  the  z  scan  peaks  from  Fig.  5  for  2.6  mW  (a)  maximum  and  (b)  minimum  and  48  mW  (c)  maximum  and  (d) 
minimum. 


radius  of  0.5  mm  at  the  cornea,  giving  an  ED60  of  approxi¬ 
mately  72  J/cm2.  The  point  in  time  at  which  each  of  the 
0.5,  0.75,  1.25,  1,75,  and  2.25  mm  beam  radii  at  the  cor¬ 
nea  begins  to  exceed  the  72  J/cm2  corneal  ED50  was  0.19, 
0.42,  1.18,  2.31,  and  3.82  s,  respectively. 

Model  results  for  short  and  long  time  intervals  are 
shown  for  1318  and  1150  nm  CW  laser  radiation  for  vari¬ 
ous  powers  delivered  to  the  cornea  in  Figs.  9-11.  Param- 


0  0.05  0.1  0.15  0.2  0.25 

Time  (s) 

Fig.  7.  Increase  in  the  beam  radius  at  the  retina  as  a  function  of 
time  for  selected  powers  delivered  to  the  cornea.  Parameters 
were  for  1318  nm  CW  laser  radiation  into  a  human  eye,  with  a 
1/e2  Gaussian  beam  radius  of  2.12  mm  at  the  virtual  cornea.  The 
y  axis  was  normalized  to  98  /im,  which  was  the  beam  radius  at 
the  virtual  retina  determined  at  time  zero. 


eters  were  for  1318  and  1150  nm  CW  laser  radiation,  in 
Figs.  9  and  11  and  in  Fig.  10,  respectively,  into  a  human 
eye  with  a  1/e2  Gaussian  beam  radius  of  2.12  mm  at  the 
virtual  cornea.  Unless  specifically  stated  in  the  caption, 
all  parameters  in  the  remaining  model  results  were  the 
same,  save  for  the  absorption  coefficients  and  index  of  re¬ 
fraction  due  to  dispersion.  In  Fig.  11,  all  parameters  were 
identical  to  those  used  to  produce  the  model  results  in 
Fig.  9,  except  the  value  of  dn/dT  was  doubled  from 
-4.46  x  10-5  to  -8.92  x  10"5  K_1. 


0  0.05  0.1  0.15  0.2  0.25 

Time  (s) 

Fig.  8.  Increase  in  the  beam  radius  at  the  retina  as  a  function  of 
time  for  selected  input  beam  radii  at  the  cornea. 
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(a)  (b) 


Fig.  9.  Increase  in  the  beam  radius  for  1318  nm  laser  radiation  at  the  retina  as  a  function  of  input  power  at  the  cornea  for  (a)  short  time 
^  1  msj  and  (b)  long  time  (1  ms  <t  =s  250  ms).  Note  the  340  W  power  at  0.350  ms  approximates  the  EDS0  threshold  level  reported  in 
Zuclich  et  al.  [34]  for  the  equivalent  wavelength.  The  y  axis  was  normalized  to  the  initial  beam  radius  of  98  pm. 


(a)  (b) 


Fig.  10.  Increase  in  the  beam  radius  for  1150  nm  laser  radiation  at  the  retina  as  a  function  of  input  power  at  the  cornea  for  (a)  short 
time  tf  =s  1  ms)  and  (b)  long  time  (1  ms  <  f  =s  250  ms).  The  y  axis  was  normalized  to  the  initial  beam  radius  of  82  pm. 


(a)  (b) 


Fig.  11.  Increase  in  the  beam  radius  for  1318  nm  laser  radiation  at  the  retina  as  a  function  of  input  power  at  the  cornea  for  (a)  short 
time  (/si  ms)  and  (b)  long  time  (1  ms  <  t  =6  250  ms)  having  doubled  the  value  of  Anl  AT  to  -8.92  X 10-6  K"1.  They  axis  was  normalized  to 
the  initial  beam  radius  of  98  pm. 
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5.  DISCUSSION 

Analysis  of  the  Schiebener  et  al.  data  for  n{T,\,P)  re¬ 
vealed  three  trends  for  water  in  the  liquid  phase  between 
283  and  363  K  at  pressures  from  0.1  to  0.2  MPa  (see  Fig. 
4  for  dn/dT  at  0.1  MPa).  First,  the  modulus  of  the 
thermo-optic  coefficient  |dn/d7’|  decreases  with  increasing 
wavelength  for  laser  radiation  between  400  and  2325  nm. 
Second,  the  modulus  of  the  thermo-optic  coefficient  in¬ 
creases  as  temperature  increases  for  liquid  water  and  fi¬ 
nally,  decreases  very  slightly  as  pressure  increases  from 
0.1  to  0.2  MPa.  This  suggests  the  magnitude  of  dn/dT  for 
ocular  media  at  body  temperature  and  standard  ocular 
pressures  should  be  approximately  two  times  greater 
than  dn/dT  values  for  water  held  in  a  cuvette  at  room 
temperature  (see  vertical  lines  in  Fig.  4).  Understanding 
these  fundamental  trends  will  lend  insight  into  predicting 
thermal  lensing  trends  in  the  eye. 

Using  data  from  Schiebener  et  al.  as  a  guide,  an  esti¬ 
mate  for  dn/dT  at  a  wavelength  of  1313  nm  between 
290-293  K  (room  temperatures)  and  pressures  near 
1  atm  (—0.1  MPa)  was  found  to  be  between  -6.96  x  10"5 
and  -7.86  x  10'5  K"1  (see  Fig.  4).  The  dn/dT  value  of 
-4.46  X  10"5  K  1  determined  to  yield  the  best  fit  to  2-scan 
data  (used  for  modeling  results  in  Figs.  5  and  6)  comes 
close  to  this  range  of  dn/dT  values.  The  error  associated 
with  the  selected  value  for  dn/dT1  is  very  difficult  to  quan¬ 
tify  given  the  sensitivity  of  the  2 -scan  data  to  the  optical 
alignment  of  the  system.  Misalignment  in  the  optical  sys¬ 
tem  for  the  2-scan  would  introduce  aberrations  not  ac¬ 
counted  for  in  the  first-order  model,  which  considers  only 
a  Gaussian  wavefront.  Further,  it  has  been  demonstrated 
that  thermal  lensing  in  water  results  in  convective  flow 
processes  introducing  nonlinear  aberrations  into  the 
steady-state  reading  [20,35-37]. 

Though  the  standard  deviation  of  the  normalized  inten¬ 
sity  for  the  experimental  2-scan  measurements  was  ±0.02 
arbitrary  units,  there  are  limitations  of  the  closed- 
aperture  2-scan.  Measures  were  taken  to  minimize  errors 
including  routine  optical  alignment  checks  before  every 
data  collection  and  careful  knife-edge  measurements  to 
characterize  the  laser  beam  profile.  However  meticulous 
these  measures  were,  they  cannot  account  for  the  uncer¬ 
tainty  with  the  beam’s  deflection  due  to  convectional  ef¬ 
fects  projected  onto  a  stationary  aperture.  Convection 
causes  small  deflections  in  the  beam’s  alignment  along 
the  x  and  y  directions,  giving  rise  to  a  nonlinear  asymmet¬ 
ric  temperature  profile,  that  are  magnified  over  the  long 
distance  to  the  aperture  in  the  far  field.  For  conditions  in¬ 
volving  higher  input  power  and  longer  exposure  times, 
such  as  the  48  mW  case  in  Figs.  5  and  6,  the  thermal  lens 
causes  the  profile  of  the  beam  to  break  its  Gaussian  dis¬ 
tribution,  forming  a  donut  shape  [20,35-37],  The  closed- 
aperture  2-scan  cannot  account  for  or  detect  conditions 
where  the  Gaussian  wavefront  is  compromised.  Alfonso 
et  al.  reported  using  a  thin  vertical  slit  in  place  of  the  iris 
aperture  shown  in  Fig.  2  to  overcome  this  problem;  how¬ 
ever,  slit  length  was  constrained  by  the  physical  dimen¬ 
sions  of  the  detector  [37].  Given  these  limitations  of  the 
single-beam-closed-aperture  2-scan,  our  future  thermal 
lensing  investigations  will  utilize  a  variation  on  the 
2-scan  to  observe  a  nonlinear  beam  profile.  Current  alter¬ 


native  investigational  methods  are  being  evaluated  and 
will  be  reported  at  a  later  date. 

The  strength  in  the  formation  of  the  nonlinear  beam 
profile  has  been  shown  to  be  directly  related  to  the  expo¬ 
sure  length,  power  of  the  radiation  source,  and  absorptiv¬ 
ity  of  the  exposed  media  [20,35-37].  For  this  reason, 
2-scan  data  were  collected  for  low  (2.6  mW)  and  high 
(48  mW)  power.  The  condition  of  a  symmetric  beam  pro¬ 
file  is  more  likely  to  be  maintained  in  the  low-power  data 
set.  Comparing  the  low-  to  high-power  2-scan  data  [Figs. 
5(a)  and  5(b),  respectively],  it  is  evident  the  amount  of 
overestimation  by  the  model  in  the  low-power  case  is 
much  less  than  in  the  high-power  case.  The  model  in  the 
low-power  transient  result  [Fig.  6(a)]  overestimated  the 
2-scan  irradiance  data  by  less  than  1%  compared  with  the 
more  than  18%  overestimated  by  the  model  in  the  high- 
power  data  set  at  1  s  [Fig.  6(c)],  The  low-power  2-scan 
data  appear  noisier  than  the  high-power  data  because  the 
irradiance  levels  were  dipping  toward  the  lower  limit  of 
the  detector.  When  the  limitations  of  the  single-beam, 
closed-aperture  2-scan  are  considered,  the  results  from 
the  first-order,  thermal-lensing,  2-scan  model  suggests 
our  model  provides  a  good  estimate  for  the  strength  of  the 
thermal  lens  effect.  The  dn/dT  value  of  -4.46  x  10"5  K_1 
found  in  the  fit  is  a  reasonable  estimate  for  trends  in  ther¬ 
mal  lensing  effects  in  the  eye.  With  the  consideration  of 
imperfections  in  the  optical  system  and  the  formation  of 
asymmetric  beam  profiles,  it  is  prudent  to  discuss  the  un¬ 
certainty  associated  with  our  model’s  results. 

Our  first-order  model’s  limitations  include  the  fact  that 
we  assume  radial  symmetry  is  maintained  throughout 
the  exposure,  neglect  heat  conduction  along  the  z  axis, 
and  treat  the  eye  as  a  simple  system  of  homogenous  slabs 
with  a  single  curved  surface  placed  on  the  cornea.  The 
creation  of  the  asymmetric  temperature  and  nonlinear 
beam  profiles  mean  our  mathematical  solution  based  on 
radially  symmetric  profiles  becomes  invalid.  The  use  of 
the  simple  eye  model  affords  a  straightforward  propaga¬ 
tion  of  the  beam  profile,  reducing  the  complexity  of  an  al¬ 
ready  difficult  heat  conduction  problem  in  the  creation  of 
the  thermal  lens.  The  first-order  thermal  lensing  results 
from  the  simple  eye  assist  us  in  approximating  the  mag¬ 
nitude  of  the  effect  the  thermal  lens  may  have  on  the  size 
of  the  laser  beam  radius  at  the  retina  in  an  eye  exposed  to 
CW  NIR  laser  radiation.  It  should  be  stressed  that  the  eye 
is  a  more  complex  optical  system  than  the  simple  model 
presented  here.  Advanced  models  of  the  eye  can  be  found 
in  Artal  et  al. ,  Tabernero  et  al. ,  and,  for  a  case  detailing 
chromatic  correction  in  the  near  infrared,  in  Fernandez 
et  al.  [38-40] . 

More  advanced  models  of  thermal  lensing  in  ocular  me¬ 
dia  to  include  heat  conduction  along  the  2 -axis  and  the 
nonlinear  beam  profile  are  currently  underway.  As  the  fi¬ 
delity  of  the  thermal  lensing  model  is  increased,  the  com¬ 
plexity  and  computational  needs,  including  resources  and 
time,  of  the  model  are  increased.  The  first-order  model  is 
computationally  inexpensive  and  fits  our  2-scan  data  rea¬ 
sonably  well,  particularly  for  the  low-power  2-scan  experi¬ 
mental  data  in  Figs.  5(a)  and  6(a).  Regardless  of  the  first- 
order  model’s  limitations,  we  find  it  suitable  for  assessing 
trends  of  thermal  lensing  in  the  eye  caused  by  near- 
infrared  laser  radiation. 
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From  model  results  shown  in  Fig.  7,  it  is  evident  the 
magnitude  of  thermal  lensing  increases  with  power.  For 
strongly  absorbed  wavelengths,  the  fluence  at  the  retina 
can  vary  significantly  with  input  power  and  time  of  the 
exposure.  The  power  levels  presented  in  Fig,  7  were  cho¬ 
sen  to  represent  the  range  of  power — up  to  the  limit  of  the 
CW  1318  nm  laser  source — used  in  the  retinal  damage 
studies  conducted  by  Zuclich  et  al.  [3,34],  Note  at  20  W 
under  the  simulated  conditions,  the  beam’s  1/e2  radius 
doubles  from  98  (im  to  196  /im  at  approximately  56  ms 
and  is  more  than  four  times  larger  than  initial  state  after 
250  ms.  Retinal  damage  threshold  studies  were  never 
successful  on  primates  for  CW  1318  nm,  and  it  has  never 
been  understood  what  the  mechanisms  were  that  limited 
damage  [1].  Model  results  from  Fig.  7  suggest  that  the 
fluence  rate  delivered  to  the  retina  would  decrease  over 
time,  potentially  preventing  damage  to  the  retina. 

The  input  beam  diameter  significantly  influences  the 
thermal  lensing  effect  as  illustrated  in  Fig.  8.  The  trend  is 
revealed  as  a  constant-power  Gaussian  input  beam  de¬ 
creases  in  radius  at  the  cornea:  the  magnitude  of  the  ther¬ 
mal  lens  greatly  increases  as  a  result  of  the  increased 
heat  source  along  the  optical  axis.  This  is  not  surprising 
given  the  physics  of  the  z-scan.  The  most  dramatic  impact 
of  the  thermal  lens  occurs  inside  the  focus  and  is  even  de¬ 
scribed  as  the  reason  for  creating  z-scan  systems  in  which 
the  sample  is  shorter  than  the  depth  of  focus  of  the  sys¬ 
tem  [12].  Thus,  as  the  depth  of  focus  (DOF)  increases,  the 
magnitude  of  the  thermal  lens  inside  the  media  should  in¬ 
crease  as  well.  The  DOF  is  given  as  2Tr\(fk/2d)2,  where  f\ 
is  the  focal  length  of  the  virtual  eye  for  the  respective 
wavelength  \,  and  d  is  the  input  diameter  of  the  laser 
beam  at  the  virtual  cornea.  Initially,  at  time  zero,  the 
smaller  input  radii  delivered  to  the  cornea  produce  the 
smallest  laser  spot  at  the  virtual  retina.  At  some  point  in 
time,  the  conditions  with  the  longer  DOF  at  the  virtual 
retina  will  give  the  larger  thermal  lens.  This  suggests 
that  when  calculations  are  made  to  increase  the  fluence 
delivered  to  the  retina  for  strongly-absorbed  laser  radia¬ 
tion,  the  condition  of  decreasing  the  input  laser  beam  ra¬ 
dius  may  in  fact  be  decreasing  the  delivered  fluence  due 
to  the  thermal  lens. 

Zuclich  et  al.  reported  for  a  1315  nm  wavelength  with  a 
350  /j.s  retinal  exposure,  the  ED50  is  0.844  J/cm2  deliv¬ 
ered  at  the  retina  [3] .  Taking  this  value  and  scaling  it  to 
the  1/e2  beam  radius  of  2.1  mm  indicates  that  0.119  J 
was  delivered  to  the  cornea.  Dividing  by  an  exposure  time 
of  350  (is  provides  an  average  power  of  340  W.  A  series  of 
powers  input  at  the  cornea  for  four  different  magnitudes, 
including  this  340  W  case,  are  shown  in  Fig.  9.  It  takes 
only  3  ms  before  the  beam  radius  at  the  virtual  retina 
doubles  in  size.  As  power  increases,  the  thermal  lens  will 
develop  faster.  Zuclich  et  al.  were  able  to  determine  reti¬ 
nal  damage  thresholds  for  short-pulsed  (<1  ms),  1318  nm 
laser  radiation  since  the  time  for  the  thermal  lens  to 
manifest  itself  is  significantly  longer  than  the  pulse  dura¬ 
tion  of  the  laser  source. 

In  Figs.  9  and  10,  the  model  was  run  for  four  different 
orders  of  magnitude  in  input  power  at  the  cornea  for  1318 
and  1150  nm  laser  radiation  wavelengths,  respectively,  to 
demonstrate  how  the  change  in  absorption  and  depth  of 
focus  affects  the  thermal  lens.  The  absorption  coefficients 


for  the  components  of  the  eye  in  Fig.  3  for  1150  and 
1318  nm  laser  radiation  are  given  in  Table  2.  Less  than 
2%  of  the  1318  nm  radiation  reaches  the  virtual  retina 
compared  with  the  near  7%  energy  transmitted  at 
1150  nm  [1].  Given  the  trends  in  literature  from  Schie- 
bener  et  al.  for  dn/dT,  the  value  of  dn/dT  should  be 
slightly  greater  for  1150  nm  compared  with  1318  nm. 
Even  with  the  value  of  dn/dT  being  the  same  for  both 
wavelengths,  these  results  can  provide  information  on  the 
behavior  of  the  thermal  lens.  Figures  9  and  10  show  that 
the  1150  nm  radiation  begins  to  invoke  the  thermal  lens 
more  strongly  for  times  less  than  250  ms,  but  at  the 
250  ms  point,  the  magnitude  of  the  thermal  lens  is  very 
close  to  the  model  results  of  the  1318  nm  radiation  case. 
Overall,  the  results  in  Fig.  9  and  10  suggest  there  is  no 
large  variation  in  the  degree  of  thermal  lensing  between 
1318  and  1150  nm  laser  radiation. 

Analysis  of  Fig.  4  suggests  that  the  value  of  the  dn/dT 
modulus  doubles  from  room  to  body  temperatures.  The 
simulations  run  for  Figs.  9  and  11  were  identical  save  for 
the  fact  that  dn/dT  was  -4.46  X 10'5  and  -8.92 
x  10  s  K~x,  respectively.  Doubling  the  magnitude  of 
dn/dT  caused  the  1/e2  beam-waist  radius  results  in  Fig. 
11(b)  to  be  approximately  33%  larger  at  250  ms  than  the 
comparable  results  in  Fig.  9(b).  These  results  demon¬ 
strate  that,  all  things  being  equal,  the  magnitude  of  the 
thermal  lens  does  increase,  but  is  not  proportional  to  the 
modulus  of  the  thermo-optic  coefficient. 

6.  CONCLUSION 

Our  mathematical  model  was  found  to  fit  experimental 
z-scan  data  reasonably  well  for  a  dn/dT  value  of  -4.46 
x  10~6  (K_1)  for  water  at  room  temperature  exposed  to 
1313  nm  CW  laser  radiation.  This  value  is  close  to  the 
predicted  dn/dT  values  based  on  analysis  of  the  Schie- 
bener  et  al.  data  (Fig.  4)  [16] .  The  trends  from  Fig.  4  sug¬ 
gest  the  value  of  dn/dT  for  water  at  body  temperature  is 
approximately  two  times  greater  than  the  dn/dT  value 
for  water  at  room  temperature.  This  conservative  value  of 
dn/dT,  -4.46  x  10-5  K_1  lends  information  on  trends  of 
thermal  lensing  in  the  human  eye  exposed  to  CW  laser  ra¬ 
diation  in  the  1150-1350  nm  region.  Our  model  provides 
a  reasonable  first-order  approximation  of  the  trends  of  the 
thermal  lensing  effect  on  the  fluence  rate  delivered  to  the 
retina.  The  thermal  lens  causes  the  fluence  rate  delivered 
to  the  retina  to  decrease  over  time,  potentially  preventing 
damage  to  the  retina  for  NIR  exposures  when  sufficient 
time  is  allowed  for  thermal  lensing  to  occur.  Modeling  re¬ 
sults  show  the  thermal  lens  takes  several  ms  to  develop. 
Retinal  damage  studies  conducted  by  Zuclich  et  al.  on 
nonhuman  primates  were  successful  for  exposure  dura¬ 
tions  less  than  a  ms,  but  could  never  explain  a  retinal  le¬ 
sion  for  exposures  greater  than  100  ms  [34].  Based  on  our 
results,  the  CW  1318  nm  laser  radiation  retinal  damage 
studies  for  exposures  greater  than  100  ms  may  have  been 
affected  by  a  decreasing  delivered  retinal  fluence  rate  pre¬ 
cluding  damage. 
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